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Abstract. We analyze the performance of a planar lens based on realistic negative 
index material in a generalized geometry. We demonstrate that the conventional su- 
perlens design (where the lens is centered between the object and the image) is not 
optimal from the resolution point-of-view, develop an analytical expression for the 
resolution limit of a generalized lens, use it to find the optimum lens configuration, 
and calculate the maximum absorption practical nearfield superlenses may have. We 
demonstrate that in contrast to the conventional superlens picture, planar imaging is 
typically accompanied by excitation of surface waves at both interfaces of the lens. 



Research on the properties of negative refractive index materials (NIMs) [1] is 
among the most rapidly developing topics in modern science that may potentially 
lead to a number of unique applications including high-performance imaging and 
lithographic systems, new types of radars, and transmission lines [2-7] . One of the 
most promising applications of NIMs involves the use of a planar slab of negative 
refraction material as an optically-perfect imaging instrument [2], known as a su- 
perlens. The physics behind the operation of a superlens and the laws governing 
the resolution limits of this unique system have instigated considerable contro- 
versy [2,8-10]. The recent analytical [11-15], numerical [16,17], and experimen- 
tal [6,16,17] results demonstrate that while a NIM-based system may outperform 
conventional (phase) lenses, its resolution logarithmically depends on material ab- 
sorption, limiting all practical applications of superlens to the near-field zone [15]. 
However, existing analytical results describing the superlens performance, are lim- 
ited to a single lens geometry, when a slab of negative-n (or negative-e) material is 
centered between an object and its image [2] and the effects of the superlens design 
on its resolution are yet to be understood. 

Addressing this fundamental question is the primary goal of the present Letter. 
In contrast to previous analytical work [12,15,16], here we consider a "generalized" 
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FIGURE 1. Schematic geometry of the generalized planar lens (a), and of the optimal planar 
lens configuration (b) as described in the text. The image is positioned at z = 2b; the possible 
lens configurations span a < b. 

variant of an imaging system (see Fig. 1), where the slab of "negative" material 
with dielectric permittivity e = — 1 + ie" and magnetic permeability fj, — — 1 + ifJ>" 
of thickness b is positioned at the distance a away from the object. Note that in 
contrast to diffraction theory presented in Ref. [18], here we are primarily interested 
in subwavelength resolution. 

We derive an analytical result for resolution of a generalized planar lens, and 
demonstrate that the resolution is maximized when a = b (the configuration re- 
cently implemented in Ref. [6]). We demonstrate that the superlens becomes im- 
practical when e" > 0.3. Further, we analyze the field distribution in the system 
and similar to what in the quasistatic limit was proved analytically [19] and indi- 
cated by other near field investigations [12,16], discover a new regime when the EM 
field has its maxima at both interfaces of NIM. We show that in contrast to most 
NIM-imaging descriptions [2,11,13-15], the field structure in generic imaging sys- 
tem is strongly influenced by this new regime, while the conventional "superlens" 
picture (with intensity minimum at the front interface) is rarely realized. Finally, 
we conclude that the optimal lens configuration suggested in this Letter minimizes 
the field intensity inside the lens and correspondingly, the total absorption in the 
imaging system. 

To better illustrate the underlying physics and simplify the analytical results 
presented in this Letter, we restrict ourselves to the case of imaging of a parallel 
slot of a thickness d <C A (with A being the free-space wavelength), extended in 
the y direction at the origin of a Cartesian system, emitting radiation with TM 
polarization (Fig. 1). The straightforward generalization of results presented here 
to the case of different shapes of object, lenses, and polarizations will be presented 
elsewhere [20]. 

The planar-lens imaging can be clearly illustrated in the wavevector space 
[11,13,15]. In this approach, the monochromatic radiation emitted by the source is 
represented as a series of waves with the same frequency u but different wavevec- 
tors k = {k x , k y , k z }. The EM field at an arbitrary point in the system can be 
calculated as a series of individual waves propagated to this point. Therefore, the 
problem of imaging an arbitrary source can be reduced to the problem of finding 
the transfer function t(x, z;k x ,u), of an individual wave with fixed (k x ,u) from 
the source (origin) to the given point in the system (x,z). For TM polarization 




FIGURE 2. (color online) Transfer functions T in (a) and r out (b) obtained from exact calcula- 
tions as described in the text (curves), and from Eqs. (2-3) (symbols); optimal and symmetric 
lens configurations with various absorptions are shown: black [solid; squares]: a = b = 0.35A, 
e" = 10 _3 ,/x" = 1CT 6 ; red [dashed; stars] a = b = 0.35A, e" = //' = 1CT 6 ; green [dash-dotted; 
triangles]: a = b/2 = 0.35A, e" = 1CT 3 ,//' = 1(T 6 ; red [dash-dot-dotted; polygons]): 
a = b/2 = 0.35A, e" = fi" = 10 -6 ; horizontal lines in (a) correspond to Eq. (9) 



considered here, it is convenient to work with the y component of the magnetic 
field: 



where a(k x ) represents the wavevector spectrum of the source [for a very thin source 
d <C A, a(k x ) ~ const]. The transfer function r is equal to Tj n inside the lens and 
Tout behind it. 

According to the properties of Fourier series [21], the component with wavevector 
k x carries the information about the x-structure of a source with a typical scale of 
2 r njk x . Therefore, the information about the fine structure of the object is being 
carried in the waves with \k x \ ^> uj/c. Since k x and k z in a plane wave are connected 
through the dispersion relation k 2 x + k\ = u 2 /c 2 , these spectral components, also 
known as evanescent waves, have imaginary k z and exponentially decay away from 
the source. The suppression of the evanescent spectrum is in fact the mechanism 
behind the resolution limit of an imaging system. 

Thus, the performance limit of a NIM-based planar lens can be related to its 
transfer function r for the evanescent part of the spectrum (\k x \ ^> lu/c) [15]. To 
calculate the transfer function, we first divide the space into three regions: before 
the lens (z < a), inside the lens (a < z < a + b), and behind the lens (z > a + b). 
We then represent the field (of a component with some fixed values of k x ,u) inside 
the first region as a sum of incident and reflected waves, the field inside the second 
region as a sum of transmitted and reflected waves, and the field in the third region 
as a transmitted wave, and use the boundary conditions to find all coefficients of 
transmission and reflection. Using this field-matching technique, described in detail 
in Ref. [15], in the limit of small absorption e" for evanescent waves 

\k x \ >ti)/cwe arrive at: 
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FIGURE 3. (color online) Intensity distributions in symmetric a = 0.35A; b = 0.7A (blue, solid) 
and optimal a — b = 0.35A (red, dashed) planar lenses with the same absorption e" = fi" = 10~ 6 . 
(a) intensity distribution along the focal line x = 0; note the intensity peaks at both lens interfaces 
(the relatively small intensity peak at the front interface of optimal system [rectangle] is shown in 
the inset); upward and downward pointing arrows show the positions of the lens in the symmetric 
and optimal configuration respectively, (b) Intensity profiles at the focal planes z = 2b of imaging 
systems; dotted black line represents the source; note that the resolution of optimal system is 
twice better than that of a symmetric structure; resolution of both systems is well-described by 
Eqs. (6,7). 
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The excellent agreement between the Eqs. (2,3) and the exact solutions of Maxwell 
equations is shown in Fig. 2. 

We now derive the resolution limit A of the generalized planar lens. Using the 
properties of Fourier analysis, the spatial size of a wavepacket at the focal point 
(x = 0, z = 2b) (see Fig. 1) can be related to its spectral width 5 through the 
"uncertainty principle" 



A • 5 = 4tt£, 



(5) 



where constant £ w 0.6 depends primarily on the geometry of a source [15,21]. 
As seen from Eq. (3), any non-zero absorption in the lens material yields an ex- 
ponential cut-off of the evanescent spectrum at some critical value of \k®\ = 5/2, 
corresponding to T out (0, 2b; k®, uj) = 1/2. This last relation can be written in the 
form of a transcendental equation for the resolution limit of a generalized planar 
lens: 
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where x = \/k® 2 c 2 /uj 2 — 1 = a/£ 2 A 2 /A 2 ~ 1- To determine the optimal configura- 
tion of the planar lens, we further simplify the Eq. (6) assuming that the system 



has subwavelength resolution A C A [22]. In this case \ ~ CV^, an d Eq. (6) 
yields: 



For the case of symmetric planar lens (b = 2a) Eqs. (6,7) are identical to the ones 
previously derived in Refs. [11,12,15,16]. We also note that Eq. (7) describes the 
resolution of a near-field "poor-man" superlens, formed by a planar slab of material 
with e = -l + ie",fjL= 1 [2,6]. 

One of the main points of this manuscript is to show that the resolution of a 
planar lens is determined not by the focal distance of the lens a, but rather by its 
thickness b, which should be minimized in order to optimize the system resolution 
(note that the imaging is possible only when b > a). Thus, for the practical 
case when the minimum separation between the object and the imaging system 
a is fixed (for example due to existence of some protective layer, etc.), the best 
possible configuration corresponds to a = b (see Fig. lb). Note that this particular 
configuration solves another problem often associated with the superlens - the local 
intensity has its maximum exactly at the focal point (as opposed to the symmetric 
lens configuration, when the intensity maximum at the back interface of NIM region 
is separated from the focus [2]), making it relatively easy to bring the optical system 
"in focus" . 

A comparison of Eq. (7) and the resolution of "conventional" near-field optics 
[21] A NF ~ 2a yields the upper limit for the absorption e^ ax of practical near-field 
superlenses when A = Anf- 



e max = 2e"*= < 2e"* ~ 0.3. (8) 



As clearly seen from this last relation, the practical applications of optical sub- 
diffraction imaging and lithography is limited to non-resonant Ag-based systems, 
while more absorbing Au-, or Al- based structures, or resonant systems [23,24] will 
have a resolution below the one achievable via conventional near-field imaging or 
almost-contact (a-separated) lithography. 

Finally, we analyze the intensity pattern and absorption in a planar lens. The 
field inside the lens region (a < z < a + b) is given by Eqs. (1,2). The intensity 
pattern at the back interface of the lens reflects a well-known effect of evanescent 
spectrum restoration [2,15]: either constant, or growing exponential (depending on 
the lens geometry) for the "weakly-evanescent" part of a spectrum uj/c < \k x \ < 
followed by exponential decay for \k x \ > k x , where k x ~ — ln(0)/6, as defined by 
Eq. (6) (see Fig. 2). 

In contrast to a widely accepted point of view [2,11], but in agreement with 
results in the quasistatic regime [19], the field at the front interface (z — a), may not 
have a sharp (exponential) minimum even for the case of sub-diffraction imaging. 
Indeed, as it is clearly seen from Eq. (2), there exists some critical wavevector 
k x ~ — ln(0)/(26), such that the exponential enhancement of the waves with k x < 



\kx\ ^ K (Fig. 2) at the front interface of a lens is possible (see Fig. 3). The 
maximum value of the field can be estimated using 

7T' c = U0,a;*2,a;)~(c ff /2r/ 6 - 1 /2, (9) 

yielding l/v / 2e" for symmetric lens configuration (b = 2a), and 1/2 for the optimal 
planar lens b = a described above. Absence of the additional strong field maximum 
and correspondingly, of the additional absorption associated with such a maximum 
further illustrates optimality of the design presented here. 

The relation of this phenomenon, which can be attributed to the excitation of 
coupled surface waves at both sides of a lens (as opposed to an excitation of a 
surface wave-"anti surface wave" pair [15]) to break-up of super-imaging and the 
onset of the diffraction limit in the system will be described in detail in our future 
work [20]. 

In conclusion, we have developed an analytical approach to the resolution of the 
generalized planar lens, used this approach to find the optimal (from the resolution 
standpoint) configuration of the lens system, and derived the maximum acceptable 
loss in the "lens" material in order to achieve a resolution gain over conventional 
near-field techniques. We also developed an analytical technique to find the field 
distribution throughout the planar imaging system, and demonstrated that there 
exists an area of resonant field excitation at the front interface of a lens due to the 
emergence of a coupled surface wave mode. 
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